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Abstract 



We present results from Brueckner-Hartree-Fock calculations for /3-stable 
neutron star matter with nucleonic and hyperonic degrees of freedom, employ- 
ing the most recent parametrizations of the baryon-baryon interaction of the 
Nijmegen group. It is found that the only strange baryons emerging in (3- 
stable matter up to total baryonic densities of 1.2 fm^^ are S~ and A. The 
corresponding equations of state are then used to compute properties of neu- 
tron stars such as masses and radii. 
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I. INTRODUCTION 



The physics of compact objects hke neutron stars offers an intriguing interplay between 
nuclear processes and astrophysical observables. Neutron stars exhibit conditions far from 
those encountered on earth; typically, expected densities p of a neutron star interior are of 
the order of 10^ or more times the 'neutron drip' density ~ 4 • 10^^ g/cm^, where nuclei 
begin to dissolve and merge together. Thus, the determination of an equation of state (EoS) 
for dense matter is essential to calculations of neutron star properties. The EoS determines 
properties such as the mass range, the mass-radius relationship, the crust thickness and the 
coohng rate. The same EoS is also crucial in calculating the energy released in a supernova 
explosion. 

At densities near to the saturation density of nuclear matter, (po = 0.16 fm~^), we expect 
the matter to be composed of mainly neutrons, protons and electrons in /3-equilibrium, since 
neutrinos have on average a mean free path larger than the radius of the neutron star. The 
equilibrium conditions can then be summarized as 

l^n = l^p + l^e, Pp = Pe, (l) 

where //j and p, refer to the chemical potential and density of particle species i, respectively. 
At the saturation density of nuclear matter, po, the electron chemical potential is of the 
order ~ 100 MeV. Once the rest mass of the muon is exceeded, it becomes energetically 
favorable for an electron at the top of the e~ Fermi surface to decay into a A Fermi sea 
of degenerate negative muons starts then to develop and, consequently, the charge balance 
needs to be changed according to Pp — Pe + P^ as well as requiring that /le — Pn- 

As the density increases, new hadronic degrees of freedom may appear in addition to 
neutrons and protons. One such degree of freedom is hyperons, baryons with a strangeness 
content. Contrary to terrestrial conditions where hyperons are unstable and decay into nu- 
cleons through the weak interaction, the equilibrium conditions in neutron stars can make 
the inverse process happen, so that the formation of hyperons becomes energetically favor- 
able. As soon as the chemical potential of the neutron becomes sufficiently large, energetic 
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neutrons can decay via weak strangeness non- conserving interactions into A hyperons leading 
to a A Fermi sea with /xa = Atrf However, one expects to appear via 

e- + n^^- + Ue, (2) 

at lower densities than the A, even though S~ is more massive. The negatively charged 
hyperons appear in the ground state of matter when their masses equal /Xe + /^n, while the 
neutral hyperon A appears when its mass equals Since the electron chemical potential in 
matter is larger than the mass difference ms- — = 81.76 MeV, S~ will appear at lower 
densities than A. For matter with hyperons as well the chemical equilibrium conditions 
become 

A^A = AtHO = A^SO = /in, 

/iS+ = /ip = yUn - fJ-e- (3) 

Hyperonic degrees of freedom have been considered by several authors, mainly within 
the framework of relativistic mean field models or parametrized effective interactions 
1^, see also Balberg et al. for a recent update. Realistic hyperon- nucleon interactions 
were employed by Schulze et al. 0, in a many-body calculation in order to study the onset 
of hyperon formation in neutron star matter. In a recent work 0, they extend their work 
to study the properties of neutron stars with hyperons, paying special attention to the role 
played by three-body nucleon forces. All these works show that hyperons appear at densities 
of the order of ~ 2po- 

In Refs. the hyperon-hyperon interaction was not included. However, it is clear 
that as soon as the S~ hyperon appears, one needs to consider the interaction between 
hyperon pairs since it will influence the single-particle energy of hyperons, hence affecting 
the equilibrium conditions from that density on and the density where other hyperons (e.g. 
the A) appear. The aim of this work is thus to present results of many-body calculations 
for /3-stable neutron star matter with hyperonic degrees of freedom, employing interactions 
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which also account for strangeness S < —1. To achieve this goal, our many-body scheme 
starts with the most recent parametrization of the free baryon-baryon potentials for the 
complete baryon octet as defined by Stoks and Rijken in Ref. ||^. This entails a microscopic 
description of matter starting from realistic nucleon-nucleon, hyperon-nucleon and hyperon- 
hyperon interactions. In a recent work |^ we have developed a formalism for microscopic 
Brueckner-type calculations of dense nuclear matter that includes all types of baryon-baryon 
interactions and allows to treat any asymmetry on the fractions of the different species 
(n,p,A,S-,S°,S+,S^ and 

Here we extend the calculations of Ref. to studies of /3-stable neutron star matter. A 
brief summary of the formalism discussed in Ref. is presented in section |I[ Our results 
are shown in section |T|. In [III A| we discuss the equation of state (EoS) and the composition 



of /^-stable matter with strangeness, using various nucleonic contributions to the EoS. Based 



on the composition of matter we discuss in section [IIIB| the possible neutron star structures. 
Our conclusions are given in section |rV|. 

II. FORMALISM 

Our many-body scheme starts with the most recent parametrization of the free baryon- 
baryon potentials for the complete baryon octet as defined by Stoks and Rijken in Ref. [§. 
This potential model, which aims at describing all interaction channels with strangeness from 
S = to S = —4, is based on SU(3) extensions of the Nijmegen potential models [|lOl for 
the 5 = and S = —1 channels, which are fitted to the available body of experimental data 
and constrain all free parameters in the model. In our discussion we employ the interaction 
version NSC97e of Ref. M, since this model, together with the model NSC97f of Ref. M, 



results in the best predictions for hypernuclear observables [jTO[. For a discussion of other 
interaction models, see Refs. P,p!Tl . 



With a given interaction model, the next step is to introduce effects from the nuclear 
medium. Here we will construct the so-called G-matrix, which takes into account short- 
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range correlations for all strangeness sectors, and solve the equations for the single-particle 
energies of the various baryons self-consistently. The G-matrix is formally given by 

{B^B2\ G{u) \B,B,) = {B^B,\ V \B,B^) + 

Y: {b,b,\v \b,b,) — ^ ^ 

x{BMGiiu)\BsB,). (4) 

Here Bi represents all possible baryons n, p, A, S^, S+, S~ and H° and their 
quantum numbers such as spin, isospin, strangeness, linear momenta and orbital momenta. 
Q is the Pauli operator which allows only intermediate states B^^Bq compatible with the 
Pauli principle, and the energy variable uj is the starting energy defined by the single- 
particle energies of the incoming external particles B^B^^. The G-matrix is solved using 
relative and centre-of-mass coordinates, see e.g., Refs. [pIpT] for computational details. The 
single-particle energies are given by 

Eb^=Tb^ + Ub^+Mb^ (5) 

where Tb^ is the kinetic energy and Mb^ the mass of baryon B^. We note that the G-matrix, 
Eq. has been solved using the standard prescription (i.e. Eb = Tb + Mb) for the 
intermediate states B^Bq. The single-particle potential Ub, is defined by 

Ub^ =Re J2 (BiB.l G{uj = Eb, + ^^J \B,B,) , (6) 

where the linear momentum of the intermediate single-particle state Bj is limited by the size 
of the Fermi surface Fj for particle species Bj. The matrix element in Eq. (|^) is properly 
antisymmetrized when the species Bi and Bj are the same. This equation is displayed in 
terms of Goldstone diagrams in Fig. [1|. Diagram (a) represents contributions from nucleons 
only as hole states, while diagram (b) has only hyperons as holes states in case we have a finite 
hyperon fraction in /5-stable neutron star matter. The external legs represent nucleons and 
hyperons. Detailed expressions for the single-particle energies and the G-matrices involved 
can be found in Ref. 0. 
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The total non-relativistic energy density, e, measured with respect to the nucleon mass, 
is obtained by adding the non-interacting leptonic contribution, ei, and the baryonic contri- 
bution, Eb, the latter being obtained from the baryon single-particle potentials 



. = ., + .. = ., + 21;/ _(^M, + — + -C/,W-M„j . (7) 

The total binding energy per baryon, is then given by 

where p is the total baryonic density. The density of a given fermion species is given by 

k'^ 

P» = ^ = ^iPT , (9) 

where Xi is the fraction of particle species i and pr = p + pi is the total density which 
includes the baryonic (p) and the leptonic one (p/). 

In order to satisfy the equations for /5-stable matter summarized in Eq. (|^), we need to 
know the chemical potentials of the particles involved. In Brueckner-Hartree-Fock (BHF) 
theory the chemical potential is taken as the single particle energy at the Fermi momentum 
of the baryon, kp , which at the lowest order reads 

PB = EB{kf^) =Mb + TB{kf^) + ^) + Ul{kf^) , (10) 

where, in the last equality, the baryon single-particle potential Ub has been split into a 
contribution, f/g, coming from the nucleonic Fermi seas (p, n) and a contribution, f/^, 
coming from the hyperonic ones (S~, A, . . .). From calculations in pure nucleonic matter 
it is well known that the nucleon chemical potential obtained from of Eq. ([T0|) differs 
considerably from the thermodynamic definition 

Therefore, for the nucleons, we replace the nucleonic contribution to the chemical potential 
in Eq. (|0D, i.e. p% = Mj, + TN^kf^) + U^{kf^), by p% = dsNN/dpN, where 
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is the nucleonic contribution to the baryonic energy density including only the interaction 
between A^A^ pairs. For the hyperons, we keep the prescription of Eq. (p!oD . As shown in 
Ref. 0, these approximations amount to ignoring the weak dependence of ,Up on the 
hyperon fractions and of f/^, Uy on the nucleon ones, and are good enough as long as the 
proton and hyperon fractions keep moderately small. Using the parabolic approximation for 
Enn one obtains [0 

fJ'pApN, 13) = I^pApn, f3 = 0) - ^/^^ ± 2/3 - /^Va^^^ j EsyniipN) , (13) 

where /3 = 1 — 2pp/ is the asymmetry parameter, with Pn = Pn + Pp- The symmetry 
energy Egym can be expressed as the difference of the energy per nucleon, S, between pure 
neutron {(3 = 1) and symmetric nuclear {(3 = 0) matter: 

1 dS 

EsyUpN) = £ipN,P= 1) - £{PN,P = 0) = 2^(/^^'/^ = 1) ' (14) 

where £ is the nucleonic contribution to the total energy per nucleon, £ = Enn/pn, and 
Pp,n{PN,P = 0) is given by 

PpApn, P = 0) = £{pN, P = 0)+ PN^^^^f£^- (15) 

The many-body approach outlined above is the lowest-order BHF method extended to the 
hyperon sector. This means also that we consider only two-body interactions. However, it is 
well-known from studies of nuclear matter and neutron star matter with nucleonic degrees 
of freedom only that three-body forces are important in order to reproduce the saturation 
properties of nuclear matter, see e.g., Ref. |T^ for the most recent approach. The effect 
of nucleon three-body forces on the properties of /5-stable matter with hyperons has been 
studied in Refs. It is found that the repulsion induced by the three-body force at high 

densities enhances substantially the hyperon population and produces a strong softening of 
the EoS. In order to include such effects, we will alternatively use for the nucleonic sector, 



the EoS of Ref. |T2| (hereafter referred to as APR98), which is obtained from a variational 



calculation using the Argonne Vig nucleon-nucleon interaction [jl3| with relativistic boost 
corrections and a fitted three-body interaction model. 

In the discussions below we will thus present two sets of results for /3-stable matter, one 
where the nucleonic contributions to the self-energy of nucleons are derived from the baryon- 
baryon potential model of Stoks and Rijken P| and one where the nucleonic contributions to 
the neutron and proton chemical potentials are calculated from the parametrization of the 
APR98 EoS discussed in Eq. (49) of Ref. Jl^. We note that replacing the nucleon-nucleon 
part of the interaction model of Ref. [§] with that from the Vis nucleon-nucleon interaction 
p!3| , does not introduce large differences at the BHF level. However, the inclusion of three- 



body forces as done in Ref. [|T2| is important. Hyperonic contributions will however all be 
calculated with the baryon-baryon interaction of Stoks and Rijken We emphasize that, 
in the present work, YN as well as YY interactions are taken into account. 



III. RESULTS 
A. Equation of state and composition of /3-stable matter 

The above models for the pure nucleonic part (NSC97e and APR98) combined with 
the hyperon contribution (NSC97e) yield the composition of /3-stable matter, up to total 
baryonic density p = 1.2 fm~^, shown in Fig. In the upper panel, results for NSC97e model 
for A^A^, yA^ and YY interactions are presented. Results combining the APR98 model for 
the nucleonic sector with the NSC97e for the YN and YY interactions are shown in the 
lower panel. In both panels, solid lines correspond to the case in which all the interactions 
A^A^, YN and YY are considered. 

As can be seen by comparing the solid lines in both panels in Fig. 0, the composition of 
/3-stable neutron star matter has a strong dependence on the model used to describe the non- 
strange sector. In both cases, is the first hyperon to appear due to its negative charge. 
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Since the APR98 EoS yields a stiffer pure nucleonic matter EoS than the corresponding one 
for NSC97e, the onset of S~ for the APR98 case occurs at a smaller density (p = 0.27 fm~'^) 
than for the NSC97e case (p = 0.34 fm~'^). In both cases, as soon as the S~ hyperon appears, 
leptons tend to disappear, totally in the APR98 case (the electron chemical potential changes 
sign at p = 1.01 fm~^, signaling the appearance of positrons), whereas in the NSC97e case 
only muons disappear. The onset of A formation takes place at higher density. Recalling 
the condition for the appearance of A, /ia = /^n = /Wp + /^e-, and that the APR98 EoS 
is stiffer due to the inclusion of three-body forces, this clearly enhances the possibility of 
creating A hyperons at lower density with this interaction model with respect to the NSC97e 
case. Indeed, the APR98 model produces A hyperons from p = 0.67 fm~^ on, whereas the 
neutron chemical potential of the NSC97e model turns out to be too small to equal the A 
one in the range of densities explored. The absence of A hyperons in the NSC97e case can 
also, in addition to a softer EoS, be retraced to a delicate balance between the nucleonic and 
hyperonic hole state contributions (and thereby to features of the baryon-baryon interaction) 
to the self-energy of the baryons considered here, see diagrams (a) and (b) in Fig. |T]. Stated 
differently, the contributions from S~, proton and neutron hole states to the A chemical 
potential are not attractive enough to lower the chemical potential of the A so that it equals 
that of the neutron. Furthermore, the increase of the chemical potential of the neutron 
with density is slowed down with the NSC97e YN interaction model since contributions 
from S~ hole states to the neutron self-energy are attractive, see e.g., Ref. for a detailed 
account of these aspects of the interaction model. We note that the isospin-dependent 
component (Lane term) of the E~ single-particle potential for the new Nijmegen interaction 
is strongly attractive, as opposed to what is found fl^ for other interactions, including 



the old Nijmegen one [|T^. This in turn implies a strong attraction for (T = 3/2) 
pairs, which is 10 times that obtained for the old Nijmegen potential at saturation density. 
These differences become more noticeable as density increases: while the H^n pairs become 
increasingly more attractive with the new Nijmegen potentials (see e.g. Fig. 5 in |Q), they 
turn out to be strongly repulsive for the old one (see e.g. Figs. 1 and 2 in 0). This 
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is why in Ref. [0 the onset density for the appearance of S~ is larger than that for free 
hyperons, whereas the reverse situation is found here (compare the S~ onset point in Fig. 
^ with what would be extracted from Fig. Within our many-body approach, no other 
hyperons appear at densities below p = 1.2 fm~'^. These results differ from present mean 
field calculations where all kinds of hyperons can appear at the densities considered 

here. Although the APR98 EoS may be viewed as the currently most realistic approach 
to the nucleonic EoS, our results have to be gauged with the uncertainty in the hyperon- 
hyperon and hyperon-nucleon interactions. Especially, if the hyperon-hyperon interactions 
tend to be more attractive, this may lead to the formation of hyperons such as the A, 
S"*", S~ and S° at lower densities. The stiffness of the nucleonic contribution, together with 
the hyperon-nucleon and hyperon-hyperon interactions play crucial roles in the appearance 
of various hyperons beyond the S~. 

In order to examine the role of the hyperon-hyperon interaction on the composition of 
/9-stable neutron star matter, the dashed lines in the lower panel of Fig. ^ show the results 
of a calculation in which only the A^A^ and YN interactions are taken into account. When 
the YY interaction is switched off, the scenario described above changes only quantitatively. 
The onset point of does not change, because S~ is the first hyperon to appear and 
therefore the YY interaction plays no role for densities below this point. We note that the 
reduction of the S~ fraction compared with the case which includes the YY interaction, 
is a consequence of neglecting the strongly attractive interaction |Q, which allows 

the energy balance {nn S^p) to be fulfilled with a smaller S~ Fermi sea. In turn, the 
reduction of the S~ fraction yields a moderate increase of the leptonic content in order to 
keep charge neutrality (in fact only muons disappear now). On the other hand, a smaller 
amount of S^'s implies less S^n pairs. Recalling that the interaction is attractive in 
this model (see e.g. Fig. 7 of Ref. 0), this means that the chemical potential of the neutrons 
becomes now less attractive. As a consequence, the A hyperons appear at a smaller density 
(p = 0.65 fm^'^) and have a larger relative fraction. 

As it has been mentioned, the composition of /5-stable matter depends strongly on the 
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model used to describe the nucleonic sector. In order to study this dependence, Fig. |^ shows 
the chemical potential of the neutron and the sum /i„ + fi^- for /5-stable matter composed 
of nucleons and free hyperons for three different NN interaction models. The solid lines 
correspond to the variational calculation denoted by APR98, which uses the Argonne Vis 
interaction and includes three-body forces and relativistic boost corrections. The dashed 
and dot-dashed lines correspond to lowest-order BHF calculations using, respectively, the 
NSC97e and the Argonne Vis potentials, the latter extracted from the results of Ref. [0. 
Dotted lines denote the A and S~ masses. In this case, the onset conditions for A and S~ 
are, respectively, /i„ = rn\ and /i„ -|- /Xe- = ?^s-- As can bee seen from the figure, the onset 
points of both hyperons are different depending on the A^A^ interaction model employed. In 
the APR98 model both the S~ and A hyperons appear at lower densities than in the lowest- 
order BHF models using the Argonne Vis or the NSC97e interactions. This is a consequence 
of the different stiffness of the EoS generated by the three A^A^ interaction models. In fact, 
the NSC97e interaction gives the softest EoS and it is not even able to produce A hyperons 
in the range of densities explored. Note that the hyperon onset points determined from 
Fig. ^ differ slightly from those observed in Fig. ^ as a consequence of the effect of the YN 
and YY interactions. Since the differences are not so large one concludes that the main 
features of the composition of matter are dominated by the pure nucleonic contribution to 
the EoS. 

At present, the APR98 EoS represents perhaps the most sophisticated many-body ap- 
proach to nuclear matter. Therefore, in what follows, we restrict our results to this A^A^ 
interaction model supplemented with the NSC97e one for the hyperonic sector. 

In Fig. ^ we show the chemical potentials in /3-stable neutron star matter for different 
baryons. We note that neither the S° nor the do appear since, as seen from the figure, 
the respective stability criteria of Eq. are not fulfilled. This is due, partly, to the fact 
that none of the S°-baryon and the S^-baryon self-energies are attractive enough. A similar 
argument applies to H*^ and S^. In the latter case the mass of the particle is ~ 1315 MeV 
and an attraction of around 200 MeV would be needed to fulfill the condition fi\ = fi^o = fin 
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at the highest density explored in this work. From the figure we see, however, that the E° 
hyperon could appear at densities close to 1.3 fm~^. 

Fig. ^ shows the EoS for four different cases: pure nucleonic matter (solid line); matter 
with nucleons and free hyperons (dotted line); matter with nucleons and hyperons inter- 
acting only with nucleons (dashed line) and, finally, matter with nucleons and hyperons 
interacting both with nucleons and hyperons (long-dashed line). Each curve corresponds to 
a different composition of /3-stable neutron star matter, obtained by requiring the equilib- 
rium conditions of Eq. (|]), with the appropriate chemical potentials for each of the four 
cases. The leptonic contribution to the EoS is also included in all the cases. As can be seen 
comparing the sohd and dotted lines the appearance of hyperons leads to a considerable 
softening of the EoS. This softening is essentially due to a decrease of the kinetic energy 
because the hyperons can be accomodated in lower momentum states and in addition have 
a large bare mass. The hyperon-nucleon interaction (dashed line) has two effects. On one 
hand, for densities up to p ~ 0.72 fm~^, the YN interaction reduces the total energy per 
baryon making the EoS even softer. On the other hand, for densities higher than p = 0.72 
fm~^, it becomes repulsive and therefore, the EoS becomes slightly stiffer than that for non- 
interacting hyperons. The contribution from the hyperon-hyperon interaction (long-dashed 
line) is always attractive producing a softening of the EoS over the whole range of densities 
explored. We note that, at high densities, the combined incorporation of the repulsive YN 
and attractive YY interactions causes the EoS to become closer to that for non-interacting 
hyperons. 



B. Structure of neutron stars 

We end this section with a discussion on neutron star properties with the above equations 
of state. 

The best determined neutron star masses are found in binary pulsars and all lie in the 



range 1.35±O.O4M0 except for the nonrelativistic pulsar PSR J1012-I-5307 of mass M 
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(2.1 ± 0.8)Mq |jT9|. Several X-ray binary masses have been measured of which the heaviest 



are Vela X-1 with M = (1.9 ± 0.2)Mq and Cygnus X-2 with M = (1.78 ± 0.2)Mq [gl 
The recent discovery of high-frequency brightness oscillations in low-mass X-ray binaries 
provides a promising new method for determining masses and radii of neutron stars, see 



Ref. |2^. The kilohertz quasi-periodic oscillations (QPO) occur in pairs and are most likely 



the orbital frequencies of accreting matter in Keplerian orbits around neutron stars of mass 
M and its beat frequency with the neutron star spin. According to Zhang et al. and 



Kaaret et al. the accretion can for a few QPO's be tracked to its innermost stable orbit. 



For slowly rotating stars the resulting mass is M ^ 2.2MQ{kllz/ uqpo) ■ For example, the 
maximum frequency of 1060 Hz upper QPO observed in 4U 1820-30 gives M ~ 2.25M0 after 
correcting for the neutron star rotation frequency. If the maximum QPO frequencies of 4U 
1608-52 {vqpo = 1125 Hz) and 4U 1636-536 {vqpo = 1228 Hz) also correspond to innermost 
stable orbits the corresponding masses are 2.IM0 and 1.9M0. These observations define a 
range of variation for the mass between M ~ 1.35M0 and M ~ 2.2Mq, which severely 
restricts the EoS for dense matter. 

In order to obtain the radius and mass of a neutron star, we have solved the Tolman- 
Oppenheimer-Volkov equation with and without rotational corrections, following the ap- 
proach of Hartle [^], see also Ref. ||14 . 



Following the discussion of Figs. H, § and ^, we use the pure nucleonic matter EoS of 



Akmal et al. [|T2[, including both relativistic boost corrections and three-body interactions. 
Relating to the discussion of the previous subsection, we study two additional CcLSGS cLS well. 
One case where we just include the hyperon-nucleon interaction as done in e.g., Refs. 
and another case where we include both the hyperon-nucleon and the hyperon-hyperon 
interactions. 

In Fig. ^ we show the resulting mass as function of these equations of state without 
rotational corrections. In Fig. we include rotational corrections to the mass. The mass- 
radius relations for the same equations of state are shown in Fig. ^. 

For /3-stable pure nucleonic matter the EoS is rather stiff compared with the EoS obtained 
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with hyperons, see Fig. |^. The EoS yields a maximum mass M = 1.89Mq without rotational 
corrections and M = 2.11Mq when rotational corrections are included. The results for the 
masses as functions of central density pc are shown by the solid lines in Figs. ^ and 0. The 
corresponding mass-radius relation for masses without rotational corrections are shown in 
Fig. g 

The second set of results refers to the case where we allow for the presence of hyperons 
and consider the hyperon-nucleon interaction but explicitly exclude the hyperon-hyperon 
one. The corresponding results are shown by the short-dashed lines in Figs. |^ and |^. 
Without rotational corrections, we obtain a maximum mass M = IAIMq whereas the 
rotational correction increases the mass to M = l.GOM©. Thus, the inclusion of the the YN 
interaction with the corresponding formation of S~ and A leads to a reduction of the mass 
by ~ 0.4 — O.SMq. This large reduction is mainly a consequence of the strong softening of 
the EoS due to appearance of hyperons. 

The last EoS employed is that which combines the nucleonic part of Ref. ||T2| with 
the computed hyperon contribution including both the hyperon-nucleon and the hyperon- 
hyperon interactions. These results are shown by the long-dashed lines in Figs. ^, |^ and |^. 
The inclusion of the hyperon-hyperon interaction leads to a further softening of the EoS in 
Fig. I, and, as can be seen from Fig. ^ this leads to an additional reduction of the total 
mass. Without rotational corrections we obtain a maximum mass M ~ 1.34M0, whilst the 
rotational correction increases the mass to M ~ 1.44M0. The size of the reduction due to 
the presence of hyperons, AM ~ 0.6 — O.TMq, and the obtained neutron star masses are 
comparable to those reported by Balberg et al. 0. 

There are other features as well to be noted from Figs. ^, ^ and The EoS with hyperons 
reaches a maximum mass at a central density pc ~ 1.3 — 1.4 fm~'^. In Fig. |^ we showed that 
the only hyperons which can appear at these densities are A and S^. If other hyperons were 
to appear at higher densities, this would most likely lead to a further softening of the EoS, 
and thereby smaller neutron star masses. 

The reader should however note that our calculation of hyperon degrees freedom is based 
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on a non-relativistic BHF approach. Although the nucleonic part extracted from Ref. ||12|| , 
including three-body forces and relativistic boost corrections, is to be considered as a bench- 
mark calculation for nucleonic degrees of freedom, relativistic effects in the hyperonic cal- 
culation could result in a stiffer EoS and thereby larger mass. However, relativistic mean 
field calculations with parameters which result in a similar composition of matter as shown 
in Fig. ^, result in similar masses as those reported in Fig. |^. In this sense, our results 
with and without hyperons may provide a lower and upper bound for the maximum mass. 
This leaves two natural options when compared to the observed neutron star masses. If the 
above heavy neutron stars prove erroneous by more detailed observations and only masses 
like those of binary pulsars are found, this may indicate that heavier neutron stars simply 
are not stable which in turn implies a soft EoS, or that a significant phase transition must 
occur already at a few times nuclear saturation densities. Our EoS with hyperons would fit 
into this case, although the mass without rotational corrections is on the lower side. Else, if 
the large masses from QPO's are confirmed, then the EoS for baryonic matter needs to be 
stiffer. This would then pose a severe problem to present hadronic models since, when the 
nucleonic part of the EoS is sufficiently stiff to support large masses, one cannot avoid the 
appearance of hyperons which, in turn, produce a softening of the EoS. 

Although we have only considered the formation of hyperons in neutron stars, transitions 
to other degrees of freedom such as quark matter, kaon condensation and pion condensation 
may or may not take place in neutron star matter. We would however like to emphasize 
that the hyperon formation mechanism is perhaps the most robust one and is likely to occur 
in the interior of a neutron star, unless the hyperon self-energies are strongly repulsive due 
to repulsive hyperon-nucleon and hyperon-hyperon interactions, a repulsion which would 
contradict present data on hypernuclei I^Gj. The EoS with hyperons yields however neutron 
star masses without rotational corrections which are even below ~ IAMq. This means 
that our EoS with hyperons needs to be stiffer, a fact which may in turn imply that more 
complicated many-body terms not included in our calculations, such as three-body forces 
between nucleons and hyperons and/or relativistic effects, are needed. 
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IV. CONCLUSIONS 



Employing the recent parametrization of the free baryon-baryon potentials for the com- 
plete baryon octet of Stoks and Rijken [Q, we have performed a microscopic many-body 
calculation of the structure of /5-stable neutron star matter including hyperonic degrees of 
freedom. The potential model employed allows for the presence of only two types of hyper- 
ons up to densities of about ten times nuclear matter saturation density. These hyperons 
are and A. The interactions for strangeness S = —1, S = —2, S = —3 and S = —4 are 
not attractive enough to allow the formation of other hyperons. 

The presence of hyperons leads however to a considerable softening of the EoS, entailing a 
corresponding reduction of the maximum mass of the neutron star. With hyperons, including 
only the hyperon-nucleon interaction, we obtain maximum masses of the order M = 1. 47 Mq 
and M = I.GOMq without and with rotational corrections, respectively. The inclusion of 
the hyperon-hyperon interaction leads to a further softening of the EoS and reduces the 
obtained masses to M = 1.34M0 and M = 1.44M0 without and with rotational corrections, 
respectively. 

The corresponding numbers with just nucleonic degrees of freedom are M = 1.89M0 and 
M = 2.IIM0, showing that the reduction in mass due to hyperonic degrees of freedom is 
~ 0.5 - O.7M0. 

Our novel result is that a further softening of the EoS is obtained when including the 
effect of the YY interaction since it is attractive over the whole density range explored. 
This is mainly due to the EE interaction which is strongly enough to develop a bound state 
0. We note that the AA attraction produced by this model is only mild, not being able to 



reproduce the experimental 2A separation energy of ABaa ~ 4 — 5 MeV ||27|| . 

Wether this additional softening is realistic or not will depend on the details of the YY 
interaction that is, unfortunatelly, not well constrained at present. New data in the S = —2 
sector, either from double-A hypernuclei or from S~-atoms, are very much awaited for. 
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FIGURES 
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FIG. 1. Goldstone diagrams for the single-particle potential Ub- a) represents the contribution 
from nucleons only as hole states while b) includes only hyperons as hole states. The wavy line 
represents the G-matrix. 
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FIG. 2. Composition of /3-stable neutron star matter. In the upper panel results for the Stoks 
and Rijken potential are presented. In the lower panel the nucleonic part of the self-energy 
of the nucleons has been replaced with the EoS of Ref. |]l2|. Solid lines in upper and lower 
panel correspond to the case in which all the interactions (nucleon-nucleon, hyperon-nucleon and 
hyperon-hyperon) are considered. Dashed lines in the lower panel correspond to the case where 
the hyperon-hyperon interaction has been switched off. 
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FIG. 3. Chemical equilibrium condition for the appearence of S~ and A hyperons for the case 
of free hyperons and three different nucleon-nucleon interactions. Dotted straight lines denote the 
rest masses of the hyperons. 
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FIG. 4. Chemical potentials in /3-stable neutron star matter as functions of the total baryonic 
density p. 
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FIG. 5. Energy per baryon in /3-stable neutron star matter as function of the total baryonic den- 
sity p, for four cases: pure nucleonic matter (solid lines); matter with nucleons and non- interacting 
hyperons (dotted line); matter with nucleons and hyperons interacting only with nucleons (dashed 
line) ; and matter with nucleons and hyperons interacting with nucleons and hyperons (long-dashed 
line). The leptonic contribution to the energy per baryon is included in all cases. 
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FIG. 6. Total mass M for various equations of state. Rotational corrections are not included. 
The solid line corresponds to the case of /3-stable matter with nucleonic degrees of freedom only, the 
short-dashed line includes also the effects of the hyperon-nucleon interaction, while the long-dashed 
line includes the three types of baryon-baryon interactions. 
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FIG. 7. Same as previous figure, but with rotational corrections. 
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FIG. 8. Mass-radius relation without rotational corrections. Notations as in the previous two 
figures. 
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